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Abstract 

In this paper, we present exact shock solutions of a coupled system of delay differential 
equations, which was introduced as a traffic-flow model called the car-following model. We 
use the Hirota method, originally developed in order to solve soliton equations. The relevant 
delay differential equations have been known to allow exact solutions expressed by elliptic 
functions with a periodic boundary conditions. In the present work, however, shock solutions 
are obtained with open boundary, representing the stationary propagation of a traffic jam. 
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1 Introduction 

Studies of traffic flow including pedestrian flow connect with a wide range of social problems such as 
efficient transport, evacuation in case of emergency and city planning. One can regard traffic flows 
as compressible fluids from a macroscopic viewpoint, or as a many-body systems of driven particles 
from a microscopic viewpoint. Accordingly, several different types of models have appeared; some 
are based on hydrodynamic equations, some are coupled differential equations, and others are 
cellular automata. 

In particular, highway-like traffic is modeled as a one-dimensional system where a number of 
particles move in a definite direction interacting with each other asymmetrically. (See FigUJ) The 
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Figure 1: An illustration of one-lane traffic. It is prohibited for particles to overtake and collide. 



one-dimensional flow presents contrasting density patterns, which change quite irregularly as the 
density of particles increases, and which finally take the form of a stable traffic jam propagating 
backwards with constant speed. In earlier work on traffic flow, car-following models were introduced 
in the form of a coupled system of DDEs: 

x n (t) = F(Ax n (t - t)) (n = l,2,...), (1) 

where x n (t) (n = 1, 2, . . .) denotes the position of the nth car at time t, Ax n (t) = x n -\(t) — x n {i) 
is the distance to the next car in front (i.e. to the (n — l)th car), and F is a given function. 
The car-following model defined by ([1]) describes a situation where each car determines its velocity 
x n (t) in terms of the distance the distance that separates it from its predecessor with a delay r, 
i.e., in terms of Ax n (t — r). The delay r is small but is never negligible, as the driver needs time 
to respond to the changing traffic situation, and hence it is one of the essential elements in traffic 
modeling. The function F, which can be determined from real traffic data, provides the optimal 
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velocity for each distance to which the drivers adjust their speed. We call it the optimal velocity 
(OV) function hereafter. 

In this paper, we present exact solutions for the car-following models with OV functions of 
exponential type (the Newell model) and tanh type (the so-called tanh model). Using the Hirota 
method, which was originally developed for solving soliton equations, we obtain shock solutions 
representing a traffic jam that propagates backwards. 

The Newell and the tanh model have been studied extensively in the past and some exact 
solutions, typically involving elliptic functions, have been obtained. However, all these solutions 
were obtained by assuming a specific traveling wave Ansatz, which reduces the original partial DDE 
to an ordinary DDE. But, as the lattice intervals of t and n are different in the original partial DDE, 
the reduced ordinary equations inevitably include two incommensurate lattice intervals, which 
makes the problem extraordinarily difficult. To solve this situation, the previous authors always 
required a strict relation between the velocity of the traveling wave and the delay (it must be equal 
to 2r). This requirement then reduces the problem to a nonlinear ordinary difference equation of 
second order (with a first order differential term) , which can be associated with addition theorems 
for the elliptic functions. It seems that up to now this was the only method available to tackle 
such nonlinear delay equations. However, in the present paper, we shall demonstrate that Hirota's 
method is very well suited to treat such equations. Exactly how this method can be applied for 
solving (fTJ) will be shown in the subsequent sections. 

The main results of this paper are reported in Section II and III. In Section II, the Newell 
model will be introduced and its bilinear form and a shock solution will be presented. Section 
III is dedicated to the tanh model, where as in Section II, a shock solution is obtained from the 
bilinear formalism. Section IV is devoted to conclusions and remarks. 



2 The Newell model 

In the epoch-making work [1] , G F Newell proposed a car-following model described by the DDE 



X n (t) = V 



1 — exp 



-±(Ax n (t-r) 



L) 



(2) 



We call it the Newell model and ([2]) the Newell equation hereafter. In [1 , particular solutions for 
@ in the case r = were presented. Also, in the subsequent work by G B Whitham [2], exact 
solutions representing solitary waves and periodic waves were found for the case r ^ 0. 

The parameters in @ should be interpreted as follows: V is the maximum allowed velocity of 
a car, 7 (the slope of the OV function at Ax n — L) corresponds to the sensitivity of the driver to 
changes in the traffic situation, and L is the minimum headway. Fig[2] shows the graph of the OV 
function. The parameters can be determined from empirical data, and in [1] the Newell model, for 
certain estimated parameter values, was shown to give a reasonable fit of the empirical data. 
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Figure 2: The optimal velocity function for the Newell model. The values of the parameters are 
as follows: V — 120, 7 = 6 and L — 5, where V indicates the maximum allowed velocity, 7 is the 
derivative of F at Ax„ = L, and L is the minimum distance between cars. 
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Firstly, eliminating uniform flows for the sake of convenience, we change the dependent variable 
from x n to y n as 

x n (t) = V t - L n + y„(t), (3) 
where the velocity Vq and headway Lq satisfy the relation 



V = V 



l-exp(-^(i -i))] (4) 



that is required of a uniform solution. Then, substituting © in ^ and using Q, the Newell 
model becomes: 

Vn{t) = (V- V ) [l - exp (-^Ay n (t -r)j] . (5) 

Considering differences of (JSJ), setting s n (t) = (j /V){y n -i(t) — y n (t)}, one obtains the following 
equation, 

— s n {t) = - exp (-s n -i(t - t)) + exp(—s n (t - r)) , (6) 

where a = 7(1 - V /V) = 7 exp (-(7 /V){L - L)). 

For instance, ([6|) is known to have the following elliptic solution [2]: 

m = , 2Q sn(»T)cn(»r)dn(»r) 

" U 8 fi(l-fc 2 sn 2 (fir)sn 2 (< ? !)+(fiT)))(l-A: 2 sn 2 (f7T)sn 2 0)' lJ 

where <j> = Q(t — 2rn), sn0 = sn(0;fc), cn^> = cn(0;fc), and dn0 = dn(</>; fc) are Jacobi's elliptic 
functions with modulus k. £1 is determined by a certain transcendental equation. 

2.1 A shock solution with parameteric freedom 

Now we apply the Hirota method [31 H] to the Newell model. In order to rewrite © as a DDE 
in algebraic form, we first set ip — exp(— s n (t)). Then, ([6]) transforms into an algebraic relation 
between ip and ifi: 

1 i^ + r) 

Following the examples of soliton equations [3j |4], we set ip = g / f , where / and g are usually 
considered to be entire functions. Then, ([5]) can be written in the following form. 

1 // / // / [if «L , q , 

a f+g+ ff ' W 

where for simplicity we use the notations, / = f(n — l,i) and / ± = f(n,t ± r) to denote n- and 
t-shifts, respectively. As a result, ((9]) can be decoupled into the following system: 

f g + f + - g + f + = Kgf - 9f) 

\lf=^f + 9 + , " ~ (10) 

where A is introduced as a coupling constant. This is the so-called bilinear form of 

Next, we assume that / = 1 + cxp(an + 2bt) and g = u + uexp(an + 2bt), where a, 6, u, and 
w are constants and the factor 2 in front of b is for later convenience. By substituting / and g 
into (flO|) . the coefficients of both exponentials give rise to the constraints on the constants. In this 
case, the number of constraints is less than that of the constants and / and g can be determined 
as follows: 

/ = 1 + exp(2b(t - ™)) 

where 6 is a free parameter. We thus obtain an exact solution of ([2|): 

a sinhQr) cosh (b(t - rn)) 
Sn (t) - log - cosh(5[t _ r(n+1)]) - i 12 ) 
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This is a shock wave with velocity U = 1/t which presents a traffic jam propagating backwards; 
and accordingly open boundary conditions are assumed (or equivalently, the period of a wave type 
solution is assumed to be infinite). It is remarkable that the shock velocity U (a macroscopic value) 
should be determined by the driver's response time r (a microscopic value). At the same time, 
this condition on the velocity means that this solution does not reduce to any previously known 
solution for which one always had U = l/2r. 

The previous solution (J7j is recovered if one takes 

f = #+(4,)M<t>), 9^o ++ (WW, (13) 
where $o(<^) = $o(<^; k) denotes the elliptic theta function with modulus k. 



3 Car-following model with the tanh OV function 

In [6l [7] , two groups of authors independently found that yet another car-following model allows 
for exact solutions, also described by elliptic functions. This model is represented by the DDE 

x n {t) = i + r,tanh ~ P ) , (14) 



where £, n, p and A are parameters. Wc call the car-following model with the OV function ([14]) 
the tanh car-following model (or simply, tanh model), hereafter. 

The tanh OV function (the RHS of (TT4"|) ) was extensively investigated in the context of the 
OV model [8] [H] ■ In particular, numerical simulation showed that this OV function reproduces the 
onset and growth of traffic jams faithfully. 

For convenience and simplicity, we shift the delay from the RHS to the LHS in (fl4[) . and 
introduce the distance variables h n (t) — (Ax n (t) — p)/2A instead of the positions x n . Then, the 
tanh car-following model is represented by the equation: 

K(t + T ) = ^J [tanh/i„_i(<) - tanh/i n (t)] . (15) 

As for the equation (|15|) , some elliptic solutions are known J5[ [7] . Typical examples can be written 
in the following forms: 

h„(t) = a sn((f>) + b, he, b he, (16) 

sn(0) + a sn^(0J + a 

where 4> = Q(t — 2rn) and Q can be regarded as a free parameter, a, b, c are constants determined 
by Q, A, n and r. 

3.1 A shock solution for the tanh model 

In terms of ip = tanhft,„, (|15j) takes the form, into an algebraic form 

'<P + = j I {l-(^ + ) 2 }(t~^. (17) 



It should be noted that, if ip is a solution of (jT7|) . — ip is also a solution. As was the case of the 
Newell model, we set ip = g/f, and thereby successfully separate ([±7]) into a system of bilinear 
equations: 

/ g + .f + -9 + .f + = Hg.f-gf) , lfi , 

\ fj = ix\[in 2 - (9 + n { } 

where A is a coupling constant. This is the bilinear form of (|15[) . 

Substituting / = 1 + exp(26t — an) and g = u + v exp(2bt — an) into (JTSJ), we find 



/ = 1 + exp(26t — an) 
2bA 

77(1 - e- 2bT ) 



g=[l- - [1 + «p(26t - an)] , ( 19 ) 
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where a is given by 



bA/r] + 1 



bA/r] - 1 + e 



-2br 



From (dHJ), we finally obtain the following exact solution of ([Til) : 



Ax n (i) = p + A log 



[ 2rjsinh(6r) 



cosh ( bt 



a \ 
2 n ) 



bA 



cosh 



(b(t 



(20) 



- 1 



This is another shock wave with velocity U — 2b/a which represents a traffic jam propagating 
backwards. It should be noticed that the velocity of (PiTj) is free, whereas the velocity of the 
solutions previously reported are all fixed at l/2r [£l[7]. 

The previous solutions (|16j) are recovered if one assume / and g are the following forms: 

/ + (0) = JA(<I>) + KM4>), = Ld x {ct>) + M§ {4>), (21) 

(which correspond to solutions of the type s n (t) = a sn(</>) + b and bj (sn(</>) + a) + c,) or 

/+ (0) = Jtf ? (4>) + K$l (0) , g+ty) = M\ (0) + M0§ (0) , (22) 

(which correspond to solutions of the type b/(sn 2 (4>) + a) + c). For the solution (p?Tj) , the constants 
J, K, L, M are determied by following quadratic relations: 



Z/ = 1 



LM 



l) 2 
"o 



J 



^) X 2 



M 



/'- 



tfa(n)i?3(n) 



ja 



(23) 



For the solution (|22|) . the constants are determied by following quartic relations: 

QA'&o^' 1 'd2^ 3 



V = 
L 2 = 
LM 
M 2 



MO) 



j 



1 ~ M " ^0 

#)tf?(fi) j2 



k diwiisi) 

r , ,2 



/' 



*3 



j9 4 



l- 



JA 

k 2 i9g(fl)#?(Q) 
5 ^Pl - 



(24) 



JA 



?9 4 



In above relations and , $j{4>) — fiji&k) denotes the elliptic theta function with modulus k and 
■dj denotes $j(0). 



4 Summary and Conclusion 

In this paper, we presented the Hirota bilinear forms of two typical car-following models, and 
obtained exact shock solutions. These are the first solutions that do not assume the velocity of 
the traveling wave to be equal to 1/2t. 

The Hirota method was originally introduced to find special solutions such as soliton solutions 
for nonlinear partial differential or difference equations. We believe that the method is effective 
for solving nonlinear equations in general and delay-differential equations in particular. 
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